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ABSTRACT 



We present the results from a detailed kinematic analysis of both ground- 
based, and Hubble Space Teles cope/ 'Faint Object Camera long-slit spectroscopy 
at sub-arcsec spatial resolution, of the narrow-line region of NGC4151. In 
agreement with previous work, the extended emission gas (R >4") is found to 
be in normal rotation in the galactic plane, a behaviour that we were able to 
trace even across the nuclear region, where the gas is strongly disturbed by 
the interaction with the radio jet, and connects smoothly with the large scale 
rotation defined by the neutral gas emission. 

The HST data, at 0".029 spatial resolution, allow us for the first time to truly 
isolate the kinematic behaviour of the individual clouds in the inner narrow-line 
region. We find that, underlying the perturbations introduced by the radio 
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ejecta, the general velocity field can still be well represented by planar rotation 
down to a radius of ~ 0".5 (30 pc), distance at which the rotation curve has its 
turnover. 

The most striking result that emerges from our analysis is that the galaxy 
potential derived fitting the rotation curve changes from a "dark halo" at the 
ENLR distances to dominated by the central mass concentration in the NLR, 
with an almost Keplerian fall-off in the \"< R < ^'interval. The observed 
velocity of the gas at 0".5 implies a mass of M ~ 10 9 M within the inner 60 pc. 
The presence of a turnover in the rotation curve indicates that this central mass 
concentration is extended. The first measured velocity point (outside the region 
saturated by the nucleus) would imply an enclosed mass of ~ 5 x 10 7 M within 
R ~ 0".15 (10 pc) which represents an upper limit to any nuclear point mass. 

Subject headings: galaxies: individual (NGC4151) - galaxies:kinematics and 
dynamics - galaxies: Seyfert 

1. Introduction 

The SABab galaxy NGC4151 hosts one of the most studied active galactic nuclei 
(AGN) in the sky, with observations from the 7-ray to the MHz range (e.g., Kriss et al. 
1995; Knop et al. 1996; Warwick et al. 1996; Ulvestad et al. 1998). The AGN spectrum 
presents a broad X-ray Fe Ka line (Yaqoob et al. 1995), optical and UV variable broad 
permitted emission lines (Crenshaw et al. 1996; Kaspi et al. 1996), as well as several 
blue-shifted narrow absorption systems (Weymann et al. 1997), and strong narrow lines 
originated in a system of clouds with up to a few kpc extension (Perez et al. 1989; Yoshida 
& Ohtani 1993). The Extended Narrow Line Region (ENLR) of NGC4151 shows line 
intensities and widths as well as kinematics consistent with quiescent gas in normal rotation 
in the galactic disk, illuminated by the central source (Penston et al. 1990; Robinson et al. 
1994). Several authors have remarked on the continuity between the kinematics of the large 
scale H I 21 cm emission and that of the ENLR (Vila-Vilaro et al. 1995; Asif et al. 1997), 
but because of the contamination by the Narrow Line Region (NLR; R < A") the ground 
based data have not yielded convincing evidence for normal rotation in the circum-nuclear 
region of NGC 4151 (Schulz 1987; Mediavilla, Arribas & Rasilla 1992). 

Narrow-band Hubble Space Telescope (HST) images (Evans et al. 1993; Boksenberg 
et al. 1995) suggested that the complex morphology of the emission clouds and filaments 
in the Narrow Line Region (NLR, R < A") is shaped not by the anisotropic character of 
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the central source's radiation but mainly by the interaction between the hot plasma of 
the radio jet and the ambient gas in the disk. In a previous paper (Winge et al. 1997, 
hereafter Paper I), we presented the initial results from HST long-slit spectroscopy of 
the NLR of NGC4151 at a spatial resolution of 0".029, isolating for the first time the 
spectra of individual clouds and demonstrating the influence of the radio ejecta in both 
the local kinematic and ionization conditions of the emission gas. Several very localized 
sub-systems of both blue and red-shifted high-velocity knots were detected, and we also 
observed off-nuclear continuum emission and marked variations on the emission line ratios 
within a few pc. Such evidence indicates that the physical conditions of the emission gas in 
individual clouds are strongly influenced by local parameters, either density fluctuations or 
shock ionization, possibly both. 

In this paper we present a detailed study of the kinematics of the gas in both the 
extended and inner NLR of NGC4151 from ground-based and HST data with high spatial 
resolution that allow us to separate the underlying velocity field of the emission gas in the 
NLR from the effects of the radio jet, and to probe its connection with the large scale 
rotation of the ENLR in the galactic disk. The most striking evidence for the interaction 
of the radio jet with the ambient gas and the presence of strong shocks, observed as 
high-velocity emission knots, localized off-nuclear continuum emission, and variations in the 
emission-line ratios in scales of a few to tens of parsec, will be discussed in a forthcoming 
paper. 

For a distance to NGC4151 of 13.3 Mpc, 0".l corresponds to a linear scale of 6.4 pc in 
the plane of the sky; a value of H = 75 kms _1 Mpc _1 is assumed throughout this paper. 

2. Observations and Data Reduction 

2.1. Ground Based data 

Long-slit spectroscopic observations of NGC4151 were obtained using the IPCS 
(Boksenberg 1972; Boksenberg & Burgess 1973) at the f/15 Cassegrain focus of the 2.5m 
Isaac Newton Telescope on March 9 to 13, 1985. The wavelength coverage includes the 
[O III] A4959,5007 and FL3 lines, at an instrumental resolution of 0.75 A (45 kms -1 ) Full 
Width Half Maximum (FWHM). The spatial resolution was 0".63. A total of 17 spectra 
were obtained at PA = 48° and 138° (parallel and perpendicular to the ENLR direction, 
respectively) at several different offsets from the nucleus. Different slit widths were used, 
from 0".22 to 0".65. The slit positions are shown in Figure [l] and the log of the observations 
is given in Table [1]. 
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The individual frames were reduced using the FIGARO image processing environment 
(Shortridge 1993), and then re-binned to a resolution of 0.24 A/pixel (60 kms -1 ). The 
continuum spectra were subtracted after correction for vignetting effects, and the spatial 
distribution of emission line fluxes, central velocity and FWHM obtained by fitting Gaussian 
functions to the line profiles using the LONGSLIT spectral analysis software (Wilkins & 
Axon 1992). 

2.2. HST FOC f/48 data 

The NLR of NGC4151 was observed using the HST Faint Object Camera (FOC) f/48 
long-slit spectrograph on July 3, 1996. The slit, 0".063 x 13". 5 in size, was positioned 
along PA = 47°. The F305LP filter was used to isolate the first order spectrum which 
covers the 3650 - 5470 A interval at a 1.58 A/pixel resolution. The spatial scale is 0".0287 
per pixel and the instrumental PSF is of the order of 0".08. The observational procedure 
consisted in obtaining first an interactive acquisition (IntAq) image in 1024 x 512 zoomed 
mode with the f/48 camera through the F220W + F275W filters to allow for an accurate 
centering of the object. While the necessary offsets were calculated, a 1247 seconds, 1024 
x 512 non-zoomed mode spectrum was taken with the slit at the IntAq position. Six 697 
seconds (non-zoomed) spectra were then obtained stepping across the NLR in 0".2 intervals, 
starting 0".6 SE of the nucleus. The three resulting SE spectra were too faint to be useful. 
The slit positions, derived as described below, are shown in Figure §, superimposed on an 
archival [O III] A5007 FOC f/96 image of the galaxy and listed in Table [|. 

The spatial relation between objects in the sky is not preserved in the uncalibrated data 
from the cameras since the raw FOC images are affected by strong geometric distortion. 
This distortion comprises two components: the external or optical, due to the off-axis 
position of the detector itself; and the internal, a combination of the distortion caused by 
the magnetic focusing of the image intensifiers and the one introduced by the spectrographic 
mirror and grating. The internal contribution is by far the most important, and it is 
strongly time and format dependent. We have used both standard IRAF| procedures for 
spectroscopic data reduction as well as specific packages developed for the FOC data (at 
stsdas.hst_calib.foc) to calibrate the data. 

Initially, all frames, including those used for subsequent calibration, were geometrically 



6 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
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corrected for the optical plus focusing induced distortions using the equally spaced grid 
of reseaux marks etched onto the first photocatode in the intensifier tube. The observed 
position of the reseaux marks were measured in the internal flat-field frames bracketing 
the observations and then compared with an equally spaced artificial grid of suitable 
size (9 x 17 reseaux marks for the 1024 x 512 mode) already corrected by the inverse 
optical distortionF]. Each individual transformation was computed fitting two dimensional 
Chebyshev polynomials of 6 th order in x and y and 5 th order in the cross-terms, and applied 
to the respective science and calibration frames. The rms uncertainties in the reseaux 
positions are 0.12 - 0.20 pixel for the 512 mode, depending on the signal to noise (S/N) 
ratio of the flat-field images. 

The remaining (mirror+grating) internal distortion along the dispersion direction was 
corrected by tracing the spectra of two stars observed in the core of the globular cluster 
47 Tuc. The stars are ~ 130 pixels apart and provide a good correction for most of the 
working area of the slit. The rms of the tracing is about 0.2 pixel. The distortion along the 
spatial direction was obtained in a similar way, tracing the brightness distribution of the 
emission lines of the planetary nebula NGC6543, with residuals of 0.10 - 0.22 pixel. The 
two corrections were combined in a single calibration file and applied simultaneously to the 
science frames. 

Since the FOC spectrograph does not contain an internal reference source, the 
NGC6543 geometrically-corrected frame was used to obtain the wavelength calibration. 
Reference wavelengths were derived from ground-based observations (Perez et al., in 
preparation), which also indicate that distortions introduced by the internal velocity field 
of the nebula are negligible at the f/48 resolution (less than 0.5 A ~ 0.3 pixel). The 
two-dimensional spectrum was collapsed along the slit, and a 6 th order Legendre polynomial 
solution found for the pixel-to-wavelength transformation, with residuals of ~ 0.04 A. The 
lines were re-identified in the original frame and a bi-dimensional wavelength calibration file 
obtained. The spatially extended emission lines in the fully calibrated frame of NGC 6543 
are measured to be within 0.13 A of their reference values. 

The procedure above can be summarized as follows: 

1. optical plus main geometric distortions are corrected using the reseaux marks and left 
a residual error of 0.12 - 0.20 pixel. 

2. residual (spectrograph) distortion in the dispersion direction is corrected tracing the 
spectra of point sources and the error is ~ 0.2 pixel. 



determined from ray-tracing models of the HST and FOC optics and available within IRAF. 
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3. residual distortion along the slit is corrected tracing the emission lines of an extended 
source and the error is 0.10 - 0.22 pixel. 

4. errors from internal velocity field of the wavelength calibrator and from the calibration 
itself are less than 0.5 and 0.13 A, respectively. 

The final errors estimated for the NGC4151 spectra, combining all the above sources 
in quadrature (including the errors from the geometric corrections applied to the calibration 
files), are 0".016 and 1.1 A rms in the spatial and dispersion directions, respectively. The 
final wavelength resolution is ~ 310 kms" 1 at [O III] A5007. 

Flux calibration was obtained using the UV standard star LDS749b. The data 
frame was geometrically corrected and wavelength calibrated as above, and the spectrum 
extracted on a 16-pixel window. This was then divided by the integration time and by the 
appropriate segment of the absolute flux table, rebinned to match the wavelength interval. 
A 6 th order spline3Q was fitted to the resulting counts per flux per second spectrum, 
averaged in 60 A intervals, generating a smooth response function. The vignetting along 
the spatial direction was corrected using the model presented in the FOC Handbook (Nota 
et al. 1996), and combined with the above response function to obtain a two-dimensional 
sensitivity calibration frame. We estimate light-losses due to the small size of the slit to be 
in the order of 20% for a point source. The relative flux of lines measured within a 0".3 - 
0".6 (10 - 20 pixels) interval is believed to be accurate at a 10 - 15% level. 

The science frames were divided by the exposure time and by the composite sensitivity 
frame. The background emission was subtracted by fitting a spline3 function along 
both spatial and dispersion directions, after masking the regions with emission lines and 
continuum. We opted to keep the fitting function order as low as possible which provides 
the best subtraction over most of the frame, even when this implied an imperfect result over 
small areas, where the background was larger and/or a more rapidly varying function of 
position. An example of the resulting 2D frames is shown in Figure |3|. The plots correspond 
to a 2". 6 x 71 A (166 pc x 4250 kms -1 ) segment of the [O III] A5007 emission line centered 
on the nuclear continuum on the PA47_1 spectrum (top) and 0".6 SW from the center of 
the slit on the PA47_3 spectrum (bottom). Multiple velocity systems and a complex cloud 
structure are seen, and some features can be easily identified in the spectra show in Figure 0, 
like the broad plume 0".37 SW of the nucleus in the top image or the double-peaked feature 
and high-velocity cloud 0".2 and 1" SW, respectively, in the 0".41 NW image. 



8 The "order" of a liner or cubic spline function in the IRAF routines refers to the number of polinomial 
pieces in the sample region. 
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To accurately determine the slit positions, we compared the [O III] A5007 luminosity 
profile of each spectrum with that derived from an archive FOC f/96 image taken with 
the F501N filter. The observed spectra were convolved with the transmission curve of the 
filter and the data collapsed along the dispersion direction. The f/96 image was rebinned 
to match the f/48 spatial resolution and the spectra light profiles compared with the sum 
of two successive lines stepping across the image at PA = 47°. The final agreement is good 
within a few percent and the final uncertainty in the position corresponds to half the size 
of the slit (or one line in the image), 0".03. 

We also retrieved from the HST Archive the NGC4151 spectra obtained on July 11, 

1995 as an engineering test of the f/48 detector. The slit was positioned at PA = 40° and the 
spectra were taken using the 1024 x 256 mode. Data reduction was complicated by the lack 
of an equivalent flat-field frame (the ones interspersed with the observations were taken in 
the 1024 x 512 zoomed mode), and contemporaneous wavelength calibration observations, 
resulting in larger uncertainties in both spatial and spectral scales. Due to the higher 
background features located in the upper and middle part of the frames, as well as the 
presence of a serious "blemish" in the detector crossing over part of the [O III] A4959,5007 
emission lines, we have chosen to optimize the reduction for the [O II] A3727 region, 
resulting in final uncertainties of 0".02 and 1.24 A in the spatial and dispersion directions, 
respectively. We note, however, that several emission and kinematic features observed in 
the 1996 [O III] data can also be identified in the 1995 frames (see Section |3|). Five data 
sets contain enough signal to be useful, and their positions, derived in a similar way to the 

1996 data, are listed in Table ^| and shown in Figure |4j. 

3. Results 

3.1. Ground-based Data 

Figure [| shows a selection of the [O III] A5007 line profiles observed along the PA48 
position at various offsets from the nucleus. In the inner 5" (the NLR), up to 3 Gaussian 
components were needed to fit the line profiles, while the extended emission was well 
represented by a single Gaussian with the instrumental resolution (~ 45 kms -1 ). The 
transition between the ENLR and the inner NLR at around 4" is particularly instructive. 
Here, the extended narrow component is joined by a second blue-shifted component of ~ 
400 kms -1 FWHM. Interior to this radius the line profiles are always broad and double. 

Figure ^| presents the result of the profile analysis for the PA48 position. It can be 
seen that the velocity structure of the extended emission (represented as stars) is similar 
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to that of a typical galactic rotation curve, indicating that this emission originates from 
gas in the galactic plane. The line width of < 45 kms" 1 is also consistent with the velocity 
dispersion of gas in the disks of normal spiral galaxies (van der Kruit & Shostak 1984). This 
component has a ratio [O IIIj/H/3 ~ 9, indicating it is photoionized by the AGN continuum. 
The general impression is that the FWHM of this primary component increases steadily 
as one approaches the nucleus while the radial velocity connects almost seamlessly onto 
the velocity field of the extended narrow component, a point we shall re-emphasize when 
we discuss the FOC f/48 results. This progression is well illustrated by the line profiles 
between 1".9 and 3". 8 SW (Figure ||c and The second component (open squares) is 
much broader, FWHM ~ 600 kms -1 , is present only in the inner 5" of the emission region, 
and is systematically blueshifted with respect to the primary. The [O III]/H/3 ratio is ~ 
7. At first sight, the physical reality of the third, narrower, component (FWHM ~ 100 
kms -1 ), which is only found very close to the nucleus, with high velocity shifts with respect 
to the main rotation curve, may be thought as questionable. However, as we shall see, it 
corresponds to structures identified in the high-spatial resolution FOC f/48 data, closely 
associated with the radio jet. 

3.2. HST FOC f/48 Data 

To study the kinematics of the gas in the HST spectra the brightest emission lines, 
[O II] A3727 and [O III] A4959,5007, were extracted in spatial windows varying from 2 to 
20 pixels along the slit, and their profiles fitted using 1 to 3 Gaussian components. The 
two [O III] lines were fitted simultaneously, constraining the physical parameters of the 
corresponding components of each line to their theoretical values (3:1 intensity ratio, 48 
A separation, and same FWHM). Figure [7] shows some examples of the fits for the 1996 
[O III] region. The spectra are identified by their distance to the center of the slit, defined 
as the line passing through the nucleus perpendicular to the slit direction, with negative 
values running towards the SW. All regions plotted are located within 1".5 from the active 
nucleus and the line profile representations vary from a single component, with very little 
broadening (second spectrum at position PA47_3) to combinations of very broad or multiple 
Gaussians. 

Figure [8] shows two regions of the [O III] A4959,5007 emission in the 1995 0".14 and 
0".52 NW (PA40_3 and PA40_5) frames, which correspond to the same kinematical features 
already remarked in Figure ^, the nuclear plume and the double peaked profile. Even with 
the presence of the blemish redward of [O III] A5007, the similarity with the nearby 1996 
profiles is evident. 
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The results of the decomposition are shown in Figures [|a to |^d for the 1996 [O III] A5007 
data, where the top panels show the FWHM of the different components superimposed over 
the brightness profiles of the radio emission (from the VLA+Merlin 5 GHz radio map of 
Pedlar et al. 1993). In these figures, the data are plotted as a function of position along 
the slit. The resulting components naturally separate into three groups: a narrow (FWHM 
< 700 kms -1 ), closest to the systemic velocity and usually the brightest component in 
the fit (filled triangles); a broad (FWHM > 700 kms -1 ) base that appears in the highest 
S/N spectra and is distributed around the radio knots (stars); and the occasional narrow 
secondary component, which can be split from the main emission by as much as ± 1000 
kms -1 (open circles). The gray region is the instrumental FWHM. The seven points 
on the PA47_1 (nuclear) position plotted as open triangles will be discussed later (see 
Sectional). 



As discussed in Paper I, there is an association between the optical and the radio 
emission, in the sense that the brightest emission-line filaments surround the radio knots, 
as expected in a scenario where the plasma of the radio jet is clearing a channel in the 
surrounding medium and enhancing the line emission along its edges by compression of 
the ambient gas (Taylor, Dyson & Axon 1992; Steffen et al. 1997a ,b ). Figures |9]a - |9|d 
present the kinematic expression of this association, with broad bases and/or high velocity 
components closely associated with the radio emission, while the "main" narrow component 
follows a more ordered pattern, strikingly resembling that of a disk rotation curve. 

The lower resolution (instrumental FWHM ~ 450 kms -1 ) and lower S/N in the 
[O II] A3727 region did not allow us to isolate as many components as for [O III] A5007. 



Nevertheless, the overall behaviour of the velocity field, shown in Figures |10| and [11] for the 
1996 and 1995 data, respectively, is the same as for the [O III] A5007 emission, with the 
narrow component tracing rotation in the disk, while the broad and secondary components 
mark the interaction of the radio plasma with the ambient gas. 

Based on our observations, we propose that the emission gas in the inner 5" of 
NGC4151 is, to the first order, produced by the central source's photoionization of the 
ambient gas located in the disk of the galaxy, and therefore, its kinematic behaviour is 
mainly planar rotation, determined by the dominant potential in the nuclear region (either 
a very concentrated but extended mass distribution or a central Massive Dark Object). 
The ionizing photons escape along a broad cone which grazes the disk, as suggested by 
Pedlar et al. (1992), and Robinson et al. (1994). From our data we can see that the 
systematic outflow previously remarked on the literature (e.g., Schulz 1990) as dominating 
the kinematics of the emission gas in the NLR of NGC4151 can be understood as an 
effect of the lack of spatial resolution of the ground based data. When the spectra of the 
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individual clouds are obtained, it is possible to decouple the gas that is in general rotation 
in the disk under the gravitational influence of the central mass concentration from that 
which is being entrained and swept along by the radio plasma, as it plunges through the 
galactic disk. This scenario has already been suggested by Vila-Vilaro et al. (1995), based 
on a higher (0.34 arcsec/pixel and ~ 1" seeing) spatial resolution ground-based spectrum 
oriented along PA = 51°, where they found that the emission within 3" from the active 
nucleus can be decomposed into one blue-shifted component located mainly SW of the 
nucleus, and a second system that they describe as appearing "to link the blue-shifted 
ENLR SW of the nucleus with the red-shifted emission to the NE. It is conceivable that it 
represents a continuation of the ENLR velocity field and hence traces the galactic rotation 
curve within the NLR." Mediavilla et al. (1992) also remarked on the existence of a gradual 
connection between the kinematic and physical properties of the galactic environment and 
the (ground-based) unresolved inner NLR. 

We also note that if the scenario of an outflow along the edges of the ionization cone 
is invoked, it is necessary for the line of sight be located outside the cone or systematic 
blue-shifted components would be projected on both SW and NE sides, which is not 
observed. A narrow cone, however, implies a ~ 25° misalignment between the outer ENLR 
and the direction of the radio jet, and therefore to a geometry that is distinct from the 
"standard" unified model. Pedlar et al. (1993) first remarked on the presence of small, 
systematic changes in the jet orientation, since the position angle of the jet in the inner 3" 
SW oscillates between 254° and 263° (74° to 83°), with an approximate mirror symmetry 
to the NE side. The abrupt change of 55° in the direction of the radio emission in the 
milliarcsec scale as seen in the VLBA radio maps of Ulvestad et al. (1998) suggests that 
effects such as precession or warping are very likely present in the collimating structure 
around the central source (Pringle 1997), but the fact that the arcsecond-scale radio jet is 
aligned within 10° indicates that these effects would be transient and their impact in the 
alignment and morphology of the NLR/ENLR/radio emission are difficult to determine. 

The presence of the variable blue-shifted optical and UV absorption lines is clear 
evidence that a gas outflow is present in the nuclear region of NGC4151. The high 
resolution GHRS spectra presented by Weymann et al. (1997) show the existence of several 
distinct systems with outflow velocities from 300 to 1600 kms -1 with respect to the nucleus, 
but the material responsible for these features is very likely to be located inside 1 pc of the 
active nucleus and the actual geometry and dynamics of the clouds is not known (Espey 
et al. 1998). Although tempting, to associate the bulk kinematics of the emission-line gas 
of the NLR on scales of several to hundreds of parsecs with the nuclear outflowing gas is 
not consistent with our data, since we observe very localized clouds with both blue- and 
red-shifted emission within the velocity range quoted for the outflows above, while the bulk 



of the gas can be well described by a planar rotation model. 



4. Modeling the Rotation Curve 

To study the gas velocity field, we concentrate on the narrowest, closer to the systemic 
velocity, and most frequently the brightest, Gaussian component (the filled triangles in 
Figures to [TT| ) as representative of the "main" NLR emission. The same argument 



was used to define the ground-based rotation curve (see Section 3T). In the ENLR, the 
general behaviour of the velocity field connects smoothly with that of the neutral gas in the 
kiloparsec scale (e.g., Vila-Vilaro et al. 1995). Here, we proceed on the assumption that 
the gas represented by the "main" component of the velocity field we isolated in the HST 
data is, to a first approximation, also participating in the general rotation of the gas in the 
galactic disk, since it also connects reasonably well with the extended emission, as can be 



seen in Figure 12 



4.1. The Model 

We have used the Bertola et al. (1991) analytic expression, which provides a simple 
parametric representation for particles (gas or stars) on circular orbits in a plane, in the 
form: 

At 

W = Vsys + p - cl)p/2 (1) 

where V sys is the systemic velocity, r is the radius in the plane of the disk and A, C , and p 
are parameters that define the amplitude and shape of the curve. If v(R, is the radial 
velocity at a position (R, \1/) in the plane of the sky, where R is the projected radial distance 
from the nucleus and ^ its corresponding position angle, we have: 

TN _ _ A R cos(* — ty ) sin i cos p i 
v mod (R, * = + Jp2 V , n2 2 . W2 2 

where 

rj = [sm 2 ($ - # ) + cos 2 i cos 2 (^ - 

where i is the inclination of the disk (i = 0° for a face-on disk) and \l/ the position angle of 
the line of nodes. 



-12- 



4.2. The Fitting 

We used a Levenberg-Marquardt non-linear least-squares algorithm to fit the above 
model. The various parameters are determined simultaneously by minimizing the residuals 
Av = (v b s — v mo d), with v mo d(R, ^; A, C Q ,p, i) and v b s (R, \E r ) being the model and observed 
radial velocities at the position (R, in the plane of the sky, respectively. 

The rotation curve expressed in Equation [| while giving a simple representation of the 
gas kinematics, has a few shortcomings from the point of view of a minimization procedure. 
The two projection angles, i and \l/ can be independently determined only if data are 
available for more than one position angle in the galaxy. In a single run of the programme, 
they are also strongly dependent on the initial values provided to the algorithm. 

The parameters A and i are strongly coupled when determining the amplitude of the 
rotation curve, so equally acceptable fits can be obtained with a large A and a small i or 
vice- versa. This situation is more acute for p ~1, and the two parameters tend to decouple 
for larger values of the exponent. The parameter C Q is determined mainly by the steepness 
of the inner part of the rotation curve, implying that the use of off-nuclear slit positions 
will tend to push its value to larger radii simply by the absence of data corresponding to 
smaller values of R. On the other hand, the exponent p is determined by the outer parts of 
the rotation curve, and its value is expected to be between 1 for a "dark halo" potential, 
and 1.5 for a system with finite mass contained within the "turn-over" radius. Evidently, 
data points for larger values of R will tend to reflect the potential generated by the mass 
distribution on larger scales. 

To test the algorithm, we used the original data (with added errors of 5-10%) from 
Bertola et al. (1991) on NGC5077 and followed the same procedure described in their 
paper. Our resulting model agrees with theirs within an rms of 8 kms -1 , with the final 
parameters differing by less than 5%. Also, applying the algorithm to an artificial rotation 
curve with spatial sampling and added errors similar to those of the NGC4151 data 
retrieves the original parameter set to better than 1%. 

4-2.1. The Ground Based data set 

To derive the projection angles of the galaxy's disk, we used the ground based data 
listed on Table [1], and excluded from the fit the points within 5" from the active nucleus 
since the spatial resolution does not allow us to separate the "main" rotating component 
from the highly disturbed gas interacting with the radio jet in the NLR. 
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Initially, the data from the three slit positions at PA = 48° were fitted with A, p, C , 
i, \I/ , and V sys as free parameters, with V sys having the same value for all data points. 
Then, to take into account any residual zero point velocity offsets between the data sets, 
each position angle was separately fitted allowing V sys to vary, but keeping the remaining 
parameters fixed to the values obtained before. The resulting systemic velocity for each 
position was subtracted from the observed values and the whole data set refitted with now 
five free parameters. Once convergence was achieved, the new model was used to obtain the 
value of V sys at each slit position, and the procedure repeated. The model was found to be 
stable at the third iteration, and the velocity offsets between data sets were smaller than 10 
kms -1 . 

Finally, the effect of the initial guess on the parameters was explored running the 
program for a wide interval of Monte Carlo search for each parameter, while the others 
were given their best fit values as the initial guesses. The resulting ranges in the best-fit 
parameters are shown in Table fj. They tend to cluster in two families characterized by 
different values in the position angle of the line of nodes (^/ )- Since the least-squares 
minimization is unable to select between them, the two average models for these families 



are also listed as Models A and B and plotted in Figure [T3|. There we see that our data 
alone is not enough to allow us to discard either of these models, so we choose the value of 
\l/ as 33°. 9, which is closer to the results quoted in the literature (\I/ = 29°. 1 ± 2.4 at R 
= 10" from H I 21 cm observations - Pedlar et al. (1992); ^ D ~ 34°- 41° at R = 2". 5 from 
Ha+[N II] measurements - Mediavilla & Arribas (1995)). We stress that the final model 
obtained from the HST data is essentially insensitive to the exact choice of value within 
this narrow range, and therefore for simplicity, in the rest of this paper we present only the 
analysis carried out using Model A. 

Other than that, the final solution is stable in all parameters: the inclination of 
the disk, i = 21°, agrees very well with other photometric and kinematic determinations 
(Simkin 1975; Bosma et al. 1977; Pedlar et al. 1992); the rotation curve derived from the 
emission gas in the ENLR is essentially flat for R > 15"(~ 1 kpc), as observed from H I 
and optical data (Pedlar et al. 1992; Robinson et al. 1994; Vila-Vilaro et al. 1995; Asif 
et al. 1997); the amplitude A is typical of the values observed for normal galaxies with 
similar absolute magnitude and Hubble type (Rubin et al. 1985). On the other hand, the 
parameter C Q ~ 440 pc, which for p — 1 is the radius at which the velocity reaches 70% of 
its maximum value, is smaller than the typical value of > 1.0 kpc, indicating a large central 
mass concentration. 



Using the Model A above to obtain the predicted rotation velocity at the positions 
sampled by the PA = 138° data, we found that the data scatters around the model without 
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much evidence of an ordered velocity field. This is not surprising for the three inner 
positions (2". 5 NE, Nucleus, 2". 5 SW) where perturbations can be induced by the radio jet 
plunging through the ambient gas. For the outermost positions, the observed points present 
much shallower gradients than expected from the model, with values smaller by 15-50 
kms -1 . The presence of turbulence on the velocity field across the outer knots has been 
already noted by other authors, and explained as effects of the expansion of the ionization 
front that produces the line emission (Asif et al. 1997) or of cloud-cloud collisions where 
the gas streaming in from the leading edge of the bar meets the one trapped in the inner 
Lindblad resonance orbits (Robinson et al. 1994; Vila-Vilaro et al. 1995). 

Therefore, in agreement with previous works, we conclude that the kinematics of the 
ionized gas in the ENLR of NGC4151 beyond a distance of about 0.5 kpc is dominated 
by the general rotation of the galactic disk, with significant perturbations present in the 
emission-line knots due to a non-circular or non-planar velocity component associated with 
gas turbulence, effects of the ionization front and/or of the presence of the galactic bar. 

Under the hypothesis that the emission in the inner 5" of the NLR of NGC4151 is 
produced by gas rotating in the disk of the galaxy, ionized mainly by the central source 
and kinematically disturbed by the interaction with the radio jet, we now assume that the 
projection of the velocity field in the plane of the sky is determined by the same geometrical 
parameters (i and \I/ ) as obtained from the ENLR rotation curve, and use the same 
procedure as above to fit the FOC f/48 data. 

Notice that we do not expect the other parameters of the fit to be the same or even 
similar to those obtained for the large scale rotation, since the high spatial resolution of 
the HST data is sampling the gas whose behaviour is governed by the very inner part of 
the gravitational potential well, while the ENLR gas reacts to the mass distribution on kpc 
scales. 



4.2.2. The 1996 FOC f/48 [0 III] data set 

As before, the full data set was first fitted with A, p, C a , and V sys as free parameters, 
while i and \l/ were kept fixed at 21° and 33°. 9, respectively, and with V sys having the same 
value for all data points, then each slit position was separately fitted with the resulting 
model but allowing V sys to vary. As a double check of the result, we also obtained a fit for 
each slit position allowing all four parameters to vary, and found that the values of V sys 
for the individual "best" and constrained fits agree well within the errors. As explained 
in Macchetto et al. (1997), the repositioning of the FOC spectrographic mirror, which 
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moves between flat-field and source exposures, can cause a shift in the zero point of the 
wavelength scale. If one of the individual data sets is shifted up or down relative to the 
others, the simultaneous fit to all data will result in a weighted value for V sys , and a larger 
value for A than what would be expected from a uniform wavelength scale. We found a ~ 
10 kms" 1 shift between the systemic velocity for position PA47_1 and both PA47_2 and 
PA47_4, obtained in the same orbit, and a ~ 50 kms -1 shift between them and PA47_3, the 
IntAq spectrum. 

The resulting systemic velocity for each slit position was then subtracted from the 
observed values and the whole data set refitted with A, p, and C Q as free parameters. 
Finally, the effect of the initial guess on the parameters was explored as for the ground 
based data, and the resulting intervals in the parameters are listed in Table |], together with 
the average model ("fit") and the one obtained using only the PA47_1 spectra ("nucleus") 



which, as discussed in Section 4.2, would be more sensitive to the actual value of C Q . These 



two models are plotted in Figure [14] as full and dotted lines, respectively. Notice that the 
points represented as open triangles in Figure |9|a had been excluded from the final fit. 
Their projected position corresponds to where the radio jet crosses the slit and therefore we 
suspect their high velocity can be due to jet entrainment. If included in the fit, their effect 
is to increase the total amplitude A of the curve, with little change in the other parameters. 
The most striking result that emerges from our analysis is that the exponent p of Eq. [I 
changes from 1 to 1.5, indicating that the potential goes from "dark halo" at the ENLR 
distances to dominated by the central mass concentration in the interval 1"< R < 4", with 
the gas rotation curve presenting an almost Keplerian fall-off in this range of radii. 



4.2.3. The 1995 and 1996 FOG f/48 [O 11] data sets 



The same procedure as for the 1996 [O III] data set was repeated using the rotation 
curves obtained from the 1995 and 1996 [O II] A3727 emission-line measurements. For the 
1995 data only the two innermost positions (PA40_2 and PA40_3) were used. The three 
data sets give similar results, shown in Table and plotted in Figures O and [16] for the 



[O II] A3727 1996 and 1995 data, respectively. The NE side of the 0".57 SE rotation curve 
is systematically blueshifted with respect to the model in the inner 1".2 from the nucleus. 
This effect can be associated with the expansion of the radio component C5 (see Figure 
2 of Paper 1), which is localized just NW of the emission sampled by the slit. The low 
S/N of the 1995 [O II] A3727 data does not allow us to completely separate the disturbed 
components, but one high velocity component can be seen at the edge of this region in 
Figure [11. 



-16- 



4.3. The HST STIS slitless data 

In March 1997, slitless CCD spectra of the NLR of NGC4151 were taken with the 
Space Telescope Imaging Spectrometer (STIS) newly installed on board HST. A description 
of the data is given in Hutchings et al. (1998). Although the [O III] A5007 spectral image 
has comparable spatial resolution to our data, it suffers from considerable confusion in the 
inner 2 - 3", where the complex velocity systems lead to an ambiguity between spatial and 
velocity information. The data are, nevertheless, worthwhile for our purpose since they 
contain velocity information extending to larger distances from the nucleus and covering 
a wider area of the NLR than our spectra. We have used the data given in Table 1 of 
the Hutchings et al. paper, and the same archival WFPC2 [O III] A5007 image to derive 
the position relative to the nucleus of the individual clouds identified in their work. We 
stress that the position ascribed by us to the clouds is uncertain by at least 2 WFPC2 
pixels (O'M). The high velocity systems detected in the STIS data were excluded and the 
remaining data within 6" from the nucleus was fitted using the same procedure as in the 
previous Sections, with V sys set to zero (the velocities in the STIS data are given relative to 
the nuclear [O III] A5007 emission). 

The resulting set of parameters is essentially identical to the one obtained for the fit 
of the FOC f/48 [O III] A5007 nuclear position alone (model [O III] A5007 Nuc. in Table |): 
A = 998 kms" 1 ; p = 1.498; C a = 0".75. The bidimensional representation of these two 
fits is presented in Figure [17|, where the STIS data and corresponding model, and the 
[O III] A5007 1996 FOC f/48 data and the Nucleus model of Table |, are shown on the 
top and bottom panels, respectively. The velocity contours are to 240 kms" 1 in steps 
of 40 kms -1 to the NE and the negative equivalents to the SW. The residuals, defined as 
Vres = V rot — V mo d e i, are negative (the model overestimates the observed velocity) on the 
left, and positive (the model underestimates the observations) on the right panels. There 
is no preferential distribution of positive and negative residuals as will be expected if the 
velocity field of the gas was partly due to a large scale radial component such as a bulk 
outflow. On the other hand, the largest residuals tend to concentrate around the position 
of the brightest radio knots, indicating that the expansion of the hot plasma introduces a 
significant non-planar component to the motion of the ambient gas. 



5. Discussion 



-17- 



5.1. The rotation curve 

Our INT data confirmed the general result from both H I and other ground based 
observations that the gas at large scales is in planar rotation, dominated by the galactic 
potential. As we approach the nucleus, it is possible to see in both the INT and HST 
data the perturbations introduced by the interaction of the radio jet with the ambient gas 
superimposed on the rotation component, which however, still dominates the bulk velocity 
field. 

The combined effect of higher spatial resolution and profile decomposition is apparent 
when we compare our data with the rotation curve of Robinson et al. (1994), obtained along 
PA = 51°. The INT curve has a smaller total amplitude in the inner 5", 180 kms -1 instead 
of ~ 350 kms -1 . Furthermore, the asymmetry remarked on previously in the ground-based 
data, where the SW peak was blue-shifted by 200 kms -1 , and the NE red-shifted by 50 
kms -1 with respect to the ENLR, disappears. The line profiles shown in Figure |5] show 
that such asymmetry is created by the mixing of the high velocity clouds plus the disk 
component. The behaviour of the line F WHM also indicates that when crossing the NLR 
boundary, the single- component profile jumps from 420 kms -1 to a little over 600 kms -1 
(see Figure 16b of Robinson et al. 1994), while our "main component" is always below 300 
kms -1 . This value is also the FWHM for the rotational component in the FOC f/48 data. 

One direct consequence of the order of magnitude improvement in spatial resolution 
provided by the FOC f/48 data is that the double-peaked, blue-shifted lines previously 
observed 3 - 4" SW of the nucleus are not evidence of a systematic bulk outflow but rather 
a complex profile created by the gas in the disk plus the disturbed component entrained 
by the radio jet. The presence of both blue- and red-shifted high velocity components in 
the HST spectrum of individual clouds reinforces the argument that the interaction of 
the jet with the ISM of the host galaxy associated with the collimated radiation from the 
central source are key factors in determining the morphology of the NLR, but not its bulk 
kinematics. Using the simple assumption that the gas in the NLR is in the same plane as 
the outer ENLR, we have been able to obtain a good parametrization for the rotation curve 
within 4" of the central source. Our results indicate that the kinematics of the NLR gas in 
the interval 30 < R < 250 pc is best represented by a thin disk in rotation around the mass 
distribution contained within the turn-over radius (~ 0.5"). The observed velocity of the 
gas at this radius, taken as purely Keplerian, would imply a mass of M ~ 10 9 M within 
the inner 60 pc. 

Although our nuclear spectrum was saturated in the inner 0".3 and the dispersion of 
the data points within the turn-over radius is rather large, the velocities measured there 
are smaller than further out, an effect that cannot be ascribed to the spatial PSF, and that 
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indicates that the 10 9 M above is not in a point mass but in an extended distribution. 
If we assume spherical symmetry for such mass distribution, our first measured velocity 
points to the NE and SW of the nucleus would imply an enclosed mass of ~ 5 x 10 8 M Q 
within R ~ 0".34 (20 pc; v obs ~ lOOkms" 1 ) or ~ 5 x 10 7 M within R ~ 0".15 (10 pc; 
Vobs~ 40 kms" 1 ), respectively This value would be an upper limit to any point mass that 
could be present there. Further observations with higher (or as high) spatial resolution but 
with better sampling in the inner region of the rotation curve and higher spectral resolution 
should be able to further constrain the nature of the mass distribution. 

The presence of a large central mass concentration was already indicated by the small 
value of the C Q parameter when fitting the kinematics the ENLR gas (see section [4.2.1 ), and 



the change of behaviour of the resulting rotation curve (from flat outside to Keplerian-like 
in the inner regions) indicates that it dominates the kinematics of the NLR. This is 
confirmed by a quick analysis of the brightness profile from the archival HST/WFPC2 
images. Assuming a mass-to-light ratio of 2, we find that the stars contribute at most 
~ 2 x 10 8 M to the mass between 0".5 and 4". Comparing with a central mass of 5 x 
10 8 - 10 9 M within a 0".5 radius and considering that the circular velocity scales on the 
square root of the mass, the presence of an extended (bulge) component would cause the 
gas in the NLR to deviate from a pure Keplerian rotation curve by < 10 - 20%. This value 
is well within the observational scatter of the data even bearing in mind the considerable 
uncertainties introduced by the effect of the central source to the brightness profile. Thus 
even though the mass inside the turnover radius is not a true point source, it is clear why 
the models yielded such a closely Keplerian value of p — 1.5 in our analysis. 

The wide-angle ionizing cone model of Pedlar et al. (1993) implies that the common 
collimation axis of the ionizing cone and the radio jets is not perpendicular to the galactic 
disk. For i ~ 21°, and a bicone opening angle of ~ 130° (see Figure 8 of Boksenberg et al. 
1995), the angle between the collimation axis and the galaxy rotation axis is ~ 25°, and the 
one between the collimation axis and the line of sight, ~ 40°. Our model for the kinematics 
of the gas in the inner NLR indicates that it is still rotating in the plane of the galaxy even 
at distances of a few tens of parsecs from the nucleus, and therefore is not directly related 
to the symmetry axis of the AGN itself. Such lack of common orientation is also indicated 
by the radio observations of Weymann et al. (1997) at subparsec scales, and by the Fe Ka 
profile presented by Yaqoob et al. (1995). If it is assumed that the Fe Ka line is produced 
in the accretion disk, the models indicate that the structure within 10 3 r g , the gravitational 
radius of the black hole, is essentially face-on. A similar situation is seen in the recent HST 
infrared observations of Centaurus A (Schreier et al. 1998), where the gas disk structure at 
~ 20 pc scales is oriented along the major axis of the bulge, an indication that its geometry 
is set by the galaxy gravitational potential rather than by the symmetry of the AGN and 
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its jet. 

5.2. Evidence for rotation in other objects 

Two other AGN, NGC 1068 and Mrk3, have been spectroscopically studied with HST 
with enough detail and spatial resolution to carry out an analysis similar to that in this 
paper. While in the first object we also observe strong but localized perturbations induced 
by the interaction of the radio jet with the ambient gas superimposed in a more general 
pattern characteristic of ordered rotation (Axon et al. 1997), in Mrk3 the gas motions 
in the region co-spatial with the radio jet are clearly dominated by the expansion of the 
cocoon of hot gas shocked and heated by the radio ejecta (Capetti et al. 1999). Such 
observations are in agreement with Nelson & Whittle (1996) results, where it was argued 
that the correlations between the stellar velocity dispersion and the [O III] profile in a large 
sample of ground-based observations of Seyferts indicate that the motions in the NLR are 
predominantly gravitational in nature, with objects with linear radio sources presenting 
broader [O III] lines. 

Evidence for an underlying rotational component is also found in ground-based studies 
of several other objects, like NGC 1365 (Hjelm & Lindblad 1996), NGC 3516 (Mulchaey et 
al. 1992; Arribas et al. 1997), NGC 2992 (Marquez et al. 1998; Allen et al. 1998), with a 
general trend for the low-ionization gas to be a better tracer of the disk component, while 
the high-ionization gas presents more deviant behaviour, usually associated with outflow. 

We remark that all these studies used the centroid of the emission lines, rather 
than the individual components and, until now, the observations did not have enough 
spatial resolution to separate the individual clouds, and show whether the often noted 
double-peaked profile is seen everywhere - implying that the outflow is actually a wind; or 
localized - as expected when the gas is entrained by the radio jet. The observations are 
also naturally biased toward the brighter emission line knots, which are more likely to be 
disturbed, either by the interaction with the radio jets, tidal effects or even the central 
source radiation. Obtaining good S/N data of the central regions of other AGN with high 
spatial resolution and even moderate spectral resolution can prove to be a useful technique 
for detecting the undisturbed gas and probing the nature of the central potential. Since 
stellar dynamical methods are rendered impotent, as the absorption lines in Seyfert spectra 
are inevitably filled in by the featureless continuum, this approach may provide the only 
means of directly determining the central object (black hole) mass in currently active nuclei. 

A potentially interesting follow-up of such studies is that, if the mass of the central 
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object M is determined using the NLR gas kinematics, the presence of continuum variability 
can, in principle, provide an estimate for another key missing piece of information on AGN 
models, the accretion rate M. Using a very simple steady-state, irradiated black-body 
approximation for the thermal structure of the accretion disk, Peterson et al. (1998) 
presented evidence of a correlation between the product (M M) and the time delay between 
different continuum wavebands in NGC7469 (Wanders et al. 1997; Collier et al. 1998), 
the only source where significant delays have been detected. If such correlation is found 
to hold for other nuclei, the two methods above would provide totally independent 
estimates for each of the quantities. However, viscous dissipation time scales are too 
long to reconcile with the observed continuum lags, which forces consideration of X-ray 
irradiated or composite models (Sincell & Krolik 1998; Collier et al. 1999). Therefore, 
with an independent measure of the central object mass, obtaining M would be dependent 
on further development of the accretion disk theory itself, and to obtain it using AGN 
continuum variability would make it necessary to carry on high sampling rate, simultaneous 
mult i- wavelength campaigns on other nearby active galaxies. 



6. Summary 

We can summarize the main results of our study of both HST and ground-based 
long-slit spectra of the inner and extended NLR of NGC4151 as follows: 

1. By decomposing the [O III] A5007 line profile in multiple Gaussian components we 
were able to trace the main kinematic component of the ENLR across the nuclear 
region, connecting smoothly the emission gas system with the large scale rotation 
defined by H I observations. 

2. Individual clouds in the NLR (R < 4") are observed to be kinematically disturbed by 
the interaction with the radio jet, but underlying these perturbations the cloud system 
is moving in a pattern compatible with disk rotation. High velocity components 
(up to ± 1000 kms -1 , relative to systemic) and broad (FWHM up to 1800 kms -1 ) 
bases are detected in the [O III] A5007 profile of the brightest clouds. Such regions 
are invariably at the edge of the radio knots, and this association, together with the 
overall morphology of the velocity field, lead us to propose that the main kinematic 
system in the inner region of NGC4151 is still rotation in the plane of the disk, 
disturbed but not defined by the interaction with the radio jet and the AGN emission. 

3. Fitting a simple expression for planar rotation to the data, we find that the ENLR 
gas (R > 4") presents a kinematic behaviour consistent with and well represented 
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by rotation in the galactic disk, with characteristics similar to other normal spiral 
systems. We obtain i = 21°, and \I/ = 34° - 43° for the inclination to the line of sight 
and position angle of the line of nodes of the disk, respectively The velocity field of 
external knots at R ~ 6" and 20" transverse to the radial direction presents evidence 
of non-planar or non-circular movements, probably associated with gas turbulence 
and streaming motions along the bar. 

4. Using the same projection angles as obtained for the ENLR, the NLR emission 
component believed to represent the continuation of the disk velocity field was 
also found to be consistent with planar rotation, although disturbed by the jet, as 
expected. However, while the velocity field of the extended ENLR gas is dominated by 
the potential of the galactic bulge, presenting a flat curve at large distances, we find 
that the behaviour of the gas in the inner NLR is best represented by a Keplerian-like 
potential, with the kinematics of the gas up to 4" dominated by the ~ 10 9 M mass 
concentration located within the turn-over radius of the rotation curve, located at ~ 
0".5. Our measurements inside the turn-over radius imply that this is an extended 
distribution rather than a point mass, and, if spherical symmetry is assumed, the 
innermost observed velocity still not affected by the spatial PSF gives a 5 x 10 7 M 
mass concentrated within a radius of about 10 pc. 
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Fig. 1. — The slit positions of the ground-based spectra of Table |l] overlayed on the 
[O III] A5007 image of the ENLR of NGC4151 from Perez et al. (1989). 

Fig. 2. — [O III] FOC f/96 image showing the position of the spectrograph slit for the 1996 
spectra. The actual length of the slit is longer than represented. North is at the top and 
East is to the left. 

Fig. 3. — Surface plots of a segment of the nuclear (PA47_1, top) and 0".41 NW (PA47_3, 
bottom) FOC f/48 1996 spectra in the [O III] A5007 region. Both regions are 2".6 by 71 A 
(166 pc by 4250 kms" 1 ). 

Fig. 4. — [O III] FOC f/96 image showing the position of the spectrograph slit for the 1995 
spectra. Orientation as in Figure [|. 

Fig. 5. — Sample of INT [O III] A5007 line profiles and Gaussian fitting. From top left to 
bottom right: (a) 0".6 NW of the nucleus; (b) 0".6 SW; (c) 1".9 SW; (d) 3".8 SW; (e) 5".3 
SW; and (f) 8". 8 SW. Notice how the narrow extended component can be traced into the 
inner 5" of the NLR. 

Fig. 6. — Result of the Gaussian decomposition procedure for the [O III] A5007 line of the 
INT PA48 spectrum. The velocities derived from the central wavelength of the individual 
components are shown in the bottom panel. The top panel presents the corresponding 
FWHM, and the solid line is the brightness profile of the narrow component. The data 
are plotted as a function of the distance to the center of the slit. For this and the other 
similar figures in the paper, error bars smaller than the point size are not plotted. (Stars) 
main narrow component; (square) broad component; (open circle) high- velocity narrow 
component. 

Fig. 7. — Gaussian fitting to a sample of [O III] A4959,5007 line profiles from the FOC f/48 
1996 spectra. The spectra are identified by their distance in arc seconds to the center of the 
corresponding slit, negative values to the SW. The number after the distance corresponds 
to how many columns in the image were coadded. The histogram is the data, the thin lines 
the individual components and underlying continuum, and the dotted line represents the 
residuals of the fit. 

Fig. 8. — Same as Figure for two positions on the 1995 FOC f/48 spectra (left), and the 
regions closest to them in the 1996 data (right) 

Fig. 9. — Result of the Gaussian decomposition procedure for the [O III] A5007 line of the 
1996 FOC f/48 spectra: (a) Nucleus, (b) 0".23 NW, (c) 0".41 NW, and (d) 0".46 NW. The 
velocities derived from the central wavelength of the individual components are shown in 
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the bottom panels. The top panels present the corresponding FWHM superposed to the 
VLA+Merlin 5 GHz brightness profile at each slit position. Both are plotted as a function 
of the distance to the center of the slit. (Filled triangle) main narrow component; (open 
circle) narrow secondary component; (stars) broad component. 

Fig. 10. — Rotation curves derived from the Gaussian decomposition of the 1996 FOC f/48 
[O II] A3727 line. From top to bottom: (a) Nucleus, (b) 0".23 NW, (c) 0".41 NW, and (d) 
0".46 NW. Symbols as in Figure 0. 



Fig. 11 — Same as in Figure |T0| for the 1995 FOC f/48 [O II] A3727 line. From top to 
bottom: (a) 0".57 SE; (b) 0".14 SE; (c) 0".14 NW;, (d) 0".43 NW; (e) 0".52 NW. Symbols 
as in Figure [| 

Fig. 12. — Comparison between the ground-based ENLR rotation curve (PA48 position; 
1st component, stars; 3rd component, open circles), with the FOC f/48 data (PA47_1; 
filled triangles) and the STIS slitless data from Hutchings et al. (1998) (open squares) 
corresponding to clouds closer to PA=48°. 

Fig. 13. — Comparison of the observed radial velocities from the INT (PA=48°) data with 
the best fitting circular rotation curve from Eq. |l|, corresponding to Models A (solid line) 
and B (dotted line) on Table |3[ The vertical lines mark the region excluded from the fit due 
to the interaction of the emission gas with the radio jet. 

Fig. 14.— Comparison between the 1996 FOC f/48 [O III] A5007 data set and the best fitting 
circular rotation curve. The dotted line corresponds to the fit to the nuclear slit position 
only, and the solid line to the fit obtained using all the four slit positions simultaneously. 
Notice that the distance scale is now relative to the nucleus. 



Fig. 15.— Same as Figure |T| for the 1996 FOC f/48 [O II] A3727 data set. The solid 
line corresponds to the fit obtained with the [O III] A5007 data set (all positions), and the 
dotted line to the fit obtained using the [O II] A3727 emission from the four slit positions 
simultaneously. 

Fig. 16.— Same as Figure |l| for the 1995 FOC f/48 [O II] A3727 data set. Here the solid line 
and dotted lines correspond to the fits obtained with the 1996 [O III] A5007 and [O II] A3727 
data sets (all positions), respectively, and the open circles represent the fit obtained using the 
1995 [O II] A3727 emission from the two innermost (0".14 NW and 0".14 SE) slit positions. 

Fig. 17. — Bidimensional representations of the rotation model for the STIS data (upper 
panels) and the [O III] A5007 Nuc. model of Table | for the FOC f/48 data (bottom). North 
is up; East to the left. The velocity contours are 0, 40, 80, 120, 160, 200, 240 kms" 1 to 
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the NE and the negative equivalents to the SW. The points are the residuals, defined as 
Vres = Vrot — Vmodei, negative on the left, and positive on the right panels. 
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Table 1. Journal of Observations - Ground Based Data 





PA 


Offset 


Exp. time 


Slit 


Seeing 


Name 


(Degrees) 


from Nucleus 


(s) 


Width 


(arcsec) 


5NE 


138 


5"NE 


1000 


0".6 


2.7 


2.5 NE 




2".5 NE 


510 


0".65 


1.7 


PA138 




Nucleus 


1050 


0".62 


2.2 


2.5 SW a 




2".5 SW 


500 


0".65 


1.7 


5SW 




5"SW 


1000 


0".6 


2.7 


10 SW 




10"SW 


1000 






20 SW 




20"SW 


1000 






3NW 


48 


3"NW 


1000 






PA48 




Nucleus 


1000 


0".22 


2.6 


3 SE 




3"SE 


1000 


0".6 


2.7 



a Spectrum obtained with a 0.6 neutral density filter. 
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Table 2. Journal of Observations - HST Data 



Name Dataset Date Exp.time(s) Format Comments 



Science Data 



NGC4151 


X2PJ010CP 


1995 Jul 11 


934 


1024x256 


PA40_1 0".57 SE 




X2PJ010DT 




934 


1024x256 


PA40_2 0".14 SE 




X2PJ010IT 




598 


1024x256 


PA40_3 0".14 NW 




X2PJ010JT 




598 


1024x256 


PA40_4 0".43 NW 




X2PJ010KT 




598 


1024x256 


PA40_5 0".52 NW 


NGC4151 


X38I0108T 


1996 Jul 03 


697 


1024x512 


PA47_1 Nucleus 




X38I0109T 




697 


1024x512 


PA47_2 0".23 NW 




X38I0102T 




1247 


1024x512 


PA47_3 0".41 NW 




X38I010AT 




697 


1024x512 


PA47_4 0".46 NW 


Calibration Frames 


NGC4151 


X2PJ010GT 


1995 Jul 11 


600 


1024x512z 


internal flat 




X38I0106T 


1996 Jul 03 


600 


1024x512 


internal flat 


NGC6543 


X3BD0102T 


1996 Sep 10 


682 


1024x512 


Planetary Neb. 




X3BD0103T 




717 


1024x256 


Planetary Neb. 




X3BD0104T 




600 


1024x512z 


internal flat 




X3BD0105T 




500 


1024x512 


internal flat 


47 Tuc. 


X34I0108T 


1996 Apr 04 


477 


1024x256z 


Globular CI. 




X34I0109T 




600 


1024x256z 


internal flat 


LDS 749b 


X3L8020BT 


1996 Dec 27 


497 


1024x512 


Flux Standard 




X3L8020AT 




800 


1024x512 


internal flat 
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Table 3. Model Fitting - Ground Based a 



Parameter 


Range 


Model A 


Model B 


A (kms -1 ) 


170 - 186 


181 


182 


P 


1.0 


1.0 


1.0 


Co n 


6.4 -7.4 


7.1 


6.6 


*o (°) 


31 - 48 


33.9 


43.7 


in 


18 - 23 


21.8 


20.9 



a PA=48°data 



Table 4. Model Fitting - HST FOC f/48 a 





1996 [O III] 




1996 [O II] 


1995 [O II] 


Parameter 


Range fit 


Nuc. 


Range fit 


Range fit 


A (kms- 1 ) 


1382 - 1412 1399 


986 


1434 - 1508 1478 


1133 - 1175 1157 


V 


1.499 1.499 


1.499 


1.499 1.499 


1.498 - 1.50 1.498 


Co (") 


1.16-1.19 1.18 


0.68 


1.04-1.10 1.07 


0.84 - 0.88 0.86 



a Using # = 33°. 9 and i= 21°. 
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